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ABSTRACT 

We report the discovery by the Large Area Telescope (LAT) onboard the Fermi Gamma-ray Space 
Telescope of high-energy 7— ray emission from the peculiar quasar PMN J0948+0022 (z — 0.5846). 
The optical spectrum of this object exhibits rather narrow H/3 (FWHM(H/3)~ 1500 km s~ 4 ), weak 
forbidden lines and is therefore classified as a narrow-line type I quasar. This class of objects is thought 
to have relatively small black hole mass and to accrete at high Eddington ratio. The radio loudness 
and variability of the compact radio core indicates the presence of a relativistic jet. Quasi simultaneous 
radio-optical- X-ray and 7-ray observations are presented. Both radio and 7-ray emission (observed 
over 5-months) are strongly variable. The simultaneous optical and X-ray data from Swift show a 
blue continuum attributed to the accretion disk and a hard X-ray spectrum attributed to the jet. 
The resulting broad band spectral energy distribution (SED) and, in particular, the 7-ray spectrum 
measured by Fermi are similar to those of more powerful FSRQ. A comparison of the radio and 7-ray 
characteristics of PMN J0948+0022 with the other blazars detected by LAT shows that this source has 
a relatively low radio and 7-ray power, with respect to other FSRQ. The physical parameters obtained 
from modelling the SED also fall at the low power end of the FSRQ parameter region discussed in 
Celotti & Ghisellini (2008). We suggest that the similarity of the SED of PMN J0948+0022 to that of 
more massive and more powerful quasars can be understood in a scenario in which the SED properties 
depend on the Eddington ratio rather than on the absolute power. 

Subject headings: quasars: individual (PMN J0948+0022) - galaxies: active - gamma rays: observa- 
tions 
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1. INTRODUCTION 

It is now widely recognised that strong radio sources 
associated with active galactic nuclei (AGN) must be 
powered by collimated relativistic energy flows (Rees 
1966). The bulk Lorentz factors (T) of these flows may 
be different in different systems and at different distances 
from the active nucleus. If a blob of plasma moving 
at relativistic speed is observed at small angles to the 
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jet axis (9 < 1/T) the observed radiation is amplified 
and the timescales shortened due to relativistic effects. 
Such systems are generically called blazars (Blandford & 
Rees 1978). Blazars have been often classified in subcat- 
egories: Flat-Spectrum Radio Quasars (FSRQ), charac- 
terized by strong and broad optical emission lines, and 
BL Lac objects, when no emission lines are apparent 
above the optical/UV continuum (e.g. Urry & Padovani 
1995). 

Remarkably, the Compton Gamma-Ray Observatory 
(CGRO) together with the first generation of Chcrcnkov 
Telescopes discovered that the SED of a number of the 
brightest blazars extend to the 7-ray range, showing two 
broad components: the first one, covering radio to soft 
X-rays, is thought to be due to the synchrotron emission 
from relativistic electrons, while the second one, cover- 
ing the hard X-/7— ray band, is generally attributed to 
inverse-Compton (IC) emission. The seed photons for 
the IC process can originate from the synchrotron radi- 
ation itself (synchrotron self-Compton, SSC; e.g. Ghis- 
ellini et al. 1985) or from an external source, like the 
accretion disk, the broad-line region or a dusty torus (ex- 
ternal Compton, EC; e.g. Dormer et al. 1992, Sikora et 
al. 1994, Blazejowski et al. 2000). 

Compiling and averaging the SED of the brightest 
blazars, Fossati et al. (1998) found an interesting trend 

- the blazar sequence - whereby for sources with low 
bolometric luminosity both components peak at high fre- 
quencies (UV/ soft X-rays for synchrotron and TeV for IC 

- High-frequency peaked BL Lacs, HBL), while, for in- 
creasing luminosities, both peaks shift to lower frequen- 
cies (Low- frequency peaked BL Lacs, LBL and FSRQ). 

Ghiscllini et al. (1998) proposed to explain the se- 
quence in terms of correlation between the random 
Lorentz factor of electrons emitting at the peaks of the 
SED (7p Ca k) and the global energy density (U) in the 
comoving frame. HBL have low U and high 7 P cak 7 while 
FSRQ have high U and low 7 poa k. The sequence can 
also be interpreted in an evolutionary frame (Bottcher & 
Dermer 2002, Cavaliere & D'Elia 2002). 

It is important to stress that the selection of objects 
with which the sequence was constructed was admittedly 
biased by the available samples, within which only a 
limited number of objects had 7-ray data (Maraschi & 
Tavecchio 2001). In fact challenges have been raised to 
the validity of the sequence (for a review, see Padovani 
2007 and references therein), which however could be 
overcome (Maraschi et al. 2008, Ghiscllini & Tavec- 
chio 2008). Today, the ongoing Fermi mission is ex- 
pected to provide a deeper and unbiased survey of the 
whole 7— ray sky, compared to that available during the 
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CGi?0/EGRET era, yielding possible surprises, as we 
will show in the present work. 

Although the origin of relativistic jets is presently still 
not understood, there is increasing evidence that the 
properties of jets are related to the properties of the ac- 
cretion flow which feeds the central black hole. In the 
case of stellar mass black holes the observed phenomenol- 
ogy points to an association of jet launching with ac- 
cretion "modes" characterized by different spectral and 
timing properties (Fender & Belloni 2004). In the extra- 
galactic domain the separation of radio sources into two 
broad classes (FRI and FRII), as well as the properties 
of their respective "beamed" representatives, BL Lacs 
and FSRQ, can be basically understood within a scenario 
based on the accretion mode: in the first class accretion 
onto the central black hole is sub-critical (in Edding- 
ton units) leading to radiatively inefficient accretion flows 
and relatively weak jets, while in the second one the ac- 
cretion rate is near critical, giving rise to bright disks and 
powerful jets (Ghisellini & Celotti 2001, Maraschi 2001, 
Maraschi & Tavecchio 2003; see, however, also Blandford 
& Levinson 1995). 

In this respect, the case of radio-loud narrow line 
Seyfert 1 (NLS1) active nuclei has received increasing at- 
tention. NLS1 arc characterized by an optical spectrum 
with narrow permitted lines FWHM(H/3) < 2000 km/s, 
the ratio between [OIIIJA5007 and H/3 smaller than 3 and 
a bump due to Fell (see, e.g., Pogge 2000 for a review). 
They exhibit also prominent soft X-ray excesses. These 
properties point to very high (near Eddington) accretion 
rates and relatively low masses (10 6 — 10 8 M Q ) (Boroson 
2002; see, however, Decarli et al. 2008, Marconi et al. 
2008). Only a small percentage of NLS1 are radio-loud 
(RL = (5 1/= 4.85GHz/S , 7/=440nm) > 10) or very radio- loud 
(RL > 100) (7% and 2.5 % respectively). Their flat radio 
spectra suggest that several of them could host relativis- 
tic jets: in fact VLBI variability indicates extremly high 
brightness temperatures and in some cases superluminal 
expansion has been observed (Komossa et al. 2006, Doi 
et al. 2006). 

Recently, Yuan et al. (2008) studied a complete sample 
of radio- loud NLS1 selected from the Sloan Digital Sky 
Survey (SDSS) sample. They find that a large fraction 
of those for which X-ray data exist show broad band 
spectra similar to those of HBL, with peaks close to the 
UV band. The study of simultaneous optical/UV/X-ray 
data of a sample of radio-loud NLS1 revealed that these 
sources often display a hard X-ray component, especially 
in bright optical and ultraviolet states, thus supporting 
the possible contribution of a relativistic jet, in some 
cases similar to FSRQ (Foschini et al. 2009). 

This class of sources thus is of extreme interest for ex- 
tending the studies of the properties of relativistic jets 
to different mass and power scales. While it is clear that 
the most radio-loud NLS1 should host relativistic jets, 
the properties of such jets at high energy are essentially 
unknown. Observations of a few selected sources in the 
TeV energy range with the Whipple and HESS Cerenkov 
telescopes were performed, but yielded only upper limits 
(Falcone et al. 2004, Aharonian et al. 2008). A de- 
tection at high energy (GeV or TeV 7-rays) is essential 
to complete the knowledge about the SED, allowing to 
constrain the inverse-Compton parameters and discuss 



analogies and differences with previously known blazars. 

Here we present the first detection, by the Fermi /L AT , 
of 7— rays from one of these radio-loud NLS1, the quasar 
PMN J0948+0022 (z = 0.5846). This makes it possi- 
ble to build the first whole SEE) from radio to 7— rays 
of a radio-loud NLS1, to constrain the inverse-Compton 
emission, and to evaluate the role of this new type of 
source in the framework of the blazar sequence and evo- 
lution. The paper is organized as follows: after a short 
identity card of the source in Sect. 2, the analysis of 
Fermi /LAT, Swift, Effelsberg and Owens Valley Radio 
Observatory data is presented in Sect. 3; Sect. 4 deals 
with the SED and the model selected to fit to the broad- 
band spectrum, while the discussion and conclusions are 
given in Sect. 5. Throughout this work, we adopted a 
ACDM cosmology from the most recent WMAP results, 
which give the following values for the cosmological pa- 
rameters: h = 0.71, fl m = 0.27, £1 A = 0.73 and with the 
Hubble-Lemaitre constant Hq = lQOh km s _1 Mpc -1 
(Komatsu et al. 2009). 

2. THE SOURCE PMN J0948+0022 

First identified in the MIT-Green Bank radio survey 
at 5 GHz (Bennett et al. 1986), this quasar is part of 
the sample drawn from the SDSS and FIRST (Faint Im- 
ages of the Radio Sky at Twenty-Centimeters) in a sys- 
tematic search for radio loud NLS1 objects (Zhou et al. 
2003, Komossa et al. 2006, Zhou et al. 2006, Yuan et 
al. 2008). The flux ratio [OIII]/H/3 = 0.1 < 3 (ac- 
cording to the most recent measurements by Zhou et al. 
2006) indicates that the Balmer lines are indeed orig- 
inated from the usual broad-line region, meaning that 
this is not a Type II AGN (Zhou et al. 2003). There is 
no obscured broad-line region, as supported by the ab- 
sence of any additional absorption in the optical (see the 
analyses by Zhou et al. 2003, 2006) and X-ray spectra 
(no additional Ah is required in the fit of Swift /XKT 
data in the present work; see Sect. 3.2). A reanalysis 
of SDSS data by Zhou et al. (2006) succeeded in sepa- 
rating the H/3 into a broad and narrow components (see 
also Rodriguez- Ardila et al. 2000), where the "broad" 
component has FWHM = 1432 ± 87 km s _1 . Therefore, 
PMN J0948+0022 has a broad and a narrow line region, 
but the FWHM of the "broad" lines is small, less than 
2000 km s _1 . 

Its strong and variable radio emission (RL > 1000), 
together with a flat and inverted radio spectrum (a r — 
-0.24, with a defined as S v cx v~ a ), make PMN J0948 + 
0022 one of the few undoubtely radio-loud narrow-line 
Seyfert 1 quasars. In addition, VLBI observations at dif- 
ferent epochs revealed high brightness temperatures and 
significant flux density variations, requiring a Doppler 
factor S > 2.5 and a viewing angle < 22° (Doi et al. 
2006). Now, a 7— ray detection is reported in the first 
list of bright 7-ray sources (significance > 10ct) detected 
by Fermi /LAT (Abdo et al. 2009a) and, specifically, 
in the List of Bright AGN Sources (LB AS, Abdo et al. 
2009b). 

3. DATA ANALYSIS 

3.1. Fermi/LAT 

PMN J0948+0022 appears in the LBAS derived from 
the first three months of the Fermi/LAT all-sky survey. 
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Fig. 1. — Fermi/LAT counts map (E > 200 McV) of the region 
centered on PKS J0948+0022 with radius 10°. The pixel size is 
0°.2. The gray scale bar is in units of LAT counts integrated in the 
5-month period. Epoch of coordinates is J2000. Nearby sources 
included in the likelihood analysis are also indicated. 

It is associated (93% confidence level) with the source 
OFGL J0948.3+0019 (Abdo et al. 2009a,b). In order 
to improve the positional accuracy for this source and 
hence increase the confidence level of the association, we 
added two more months of data to the three months of 
the LB AS, thus providing 5 months of data (from August 
through December 2008). 

The data from the Large Area Telescope (LAT, At- 
wood ct al. 2008) were analyzed using the Fermi/LAT 
Science Tools v 9 . 8 . 2, a specific LAT software pack- 
age 52 . In summary, we performed the analysis using 
the guidelines described in detail in the LBAS paper 
(Abdo ct al. 2009a, b), but with some specific features 
to best fit the characteristics of this source. Events 
of "Diffuse" class in the Science Tools, coming from 
zenith angles < 105° (to avoid Earth's albedo) were ex- 
tracted from a region with a 10° radius centered on the 
coordinates of the radio position of PMN J0948+0022 
{RA = 09 h 48 m 57 s .6 and Dec = +00°22'12".0, J2000). 
Because of calibration uncertainties at low energies, data 
were selected with energies above 200 MeV, with no 
7— ray detected at energies > 3 GcV. The 7— ray back- 
ground is mainly due to three components: the diffuse 
emission from the Milky Way, the diffuse extragalactic 
background and the instrumental background. All these 
components were modeled and accounted for in the anal- 
ysis. 

An unbinncd likelihood algorithm, implemented in the 
LAT Science Tools as the gtlike task, was used to 
analyze the data. PMN J0948+0022 was modelled to- 
gether with four other nearby sources: PKS 0907+022, 
PKS 0907-023, PMN J1016+0512, PMN J0953-0840, 
to take possible contaminations into account, due to 
the strongly energy-dependent point-spread function (see 
Figure 1). 

The best fit position of the LAT source, after 5 
months of integration, is RA = 09 h 48 m 58 s and Dec = 

52 http://fermi.gsfc.nasa.gov/ssc/data/analysis/software/ 



+00°22'48" (J2000), with a 95% error radius of 12'36" 
(0°.21) and detection test statistics TS = 305 (a 17cr 
detection; see Mattox et al. 1996 for the definition of 
TS). Figure 1 shows the LAT counts map with the 
peak of counts consistent with the radio position of 
PMN J0948+0022. The Figure-of-Merit (FoM) method 
to associate 7— ray sources with known radio counter- 
parts was used (Sowards-Emmcrd ct al. 2003; Abdo et 
al. 2009a,b). The value of the FoM increased from 18.64 
and 93% confidence level (Bright AGN List, 3 months 
data) to 53.92 and 99% probability that the association 
is correct (5 months data). 

Fig. 2 (A) displays the light curve for E > 200 MeV 
integrated over 10-day time bins between 1 August 2008, 
00 UT, to 1 January 2009, 00 UT (MJD 54679 - 54832). 
The corresponding photon indices T, defined as F(E) oc 
E~ T , arc shown in Fig. 2 (B). The light curve indicates 
clear variability of a factor of <~ 10 over a timescale of 
weeks, confirmed with a x 2 — 34.5 for 14 degrees of free- 
dom (% 2 = 2.5) for a fit with a constant, and an excess 
variance of 0.19 ± 0.06 (see Nandra et al. 1997 for a defi- 
nition of excess variance). On the other hand, the photon 
index does not show evident variability and can be fitted 
with a constant (\ 2 — 0.80). The radio data at 15 GHz 
in Fig. 2 (C), described in Section 3.3, displays no vari- 
ability during the same time interval. We note, however, 
that radio data are missing during the episode of vari- 
ability at 7-rays at the beginning of the LAT lightcurve 
(Fig. 2 (A)) and a strong radio outburst was observed 
well before the launch of Fermi. 

The spectrum averaged over the whole data set was 
initially fitted with a single power-law model with T = 
2.6 ± 0.1 and flux (E > 200 MeV) equal to (4.0 ± 0.3) x 
10~ 8 ph cm" 2 s" 1 (TS = 298). A fit with a broken- 
power law model gives a slight increase of the TS. The 
spectral parameters are the following: the break energy 
is Efr = 1.0 ± 0.4 GeV, while the photon index for E < 
E h is Ti = 2.3 ± 0.2 and for E > E h is T 2 = 3.4 ± 
0.5. The integrated flux (E > 200 MeV) is (3.9 ± 0.3) x 
10~ 8 ph cm" 2 s" 1 (TS = 305). The likelihood test ratio 
gives 97% probability in favor of the broken power-law 
model with respect to the single power-law. 

It is worth noting that the quoted errors are statistical 
only. Systematic errors should be added. According to 
the studies on the Vela Pulsar (Abdo et al. 2009c), our 
current conservative estimates of systematic errors are < 
30% for flux measurements and 0.1 for the photon index. 
Significant reduction of such systematic uncertainties is 
expected once the calibration of the LAT instrument is 
completed. 

3.2. Swift 

On 5 December 2008 at 02:25 UTC (MJD 54805.10), 
Swift observed PMN J0948 + 0022 (ObsID 00031306001, 
exposure ~ 4 ks). For the screening, reduction and 
analysis of the data from the three instruments (BAT, 
Barthclmy ct al. 2005; XRT, Burrows et al. 2005; 
UVOT, Roming et al. 2005) onboard the Swift satellite 
(Gehrels et al. 2004), we used the HEASoft v. 6.6.1 
software package, together with the CALDB updated on 
9 December 2008. 

The X-Ray Telescope XRT (0.2-10 kcV energy band) 
was used in photon counting mode and no evidence of 
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Fig. 2.— (panel A) Fermi/LAT flux for E > 200 MeV in units of [10 -7 ph cm 2 s x ] with 10-day time bins and (panel B) the 
corresponding photon index, (panel C) OVRO radio flux at 15 GHz [Jy]. Time scale is from 4 September 2007 at 16:25:41 UTC to 11 
February 2009 06:53:48 UTC and the day in the X-axis refers to the beginning of LAT observations (1 August 2008). 



TABLE 1 

Summary of results from analysis of the Swift data obtained on 5 December 2008 (ObsID 00031306001). See the text for 

DETAILS. 



BAT (20 - 100 koV) 




Exposure 
[ks] 




Flux 2 0-100keV 

[10- 10 erg cm" 2 s" 1 ] 






4.3 




< 7.0 




XRT (0.2 - 10 koV) 


Exposure 

M 


[10 20 cm- 2 ] 


r 


Normalization at 1 kcV FluX2-iokcV 
[10~ 3 ph cm -2 s" 1 keV -1 ] [10~ 12 erg cm -2 s _1 ] 


X 2 /dof 


4.2 


5.22 


1.83 ± 0.17 


1.5 ±0.3 2.2 ±0.2 


5.1/8 


UVOT (observed magnitudes) 


V 


b 


u 


uvwl uvm2 


uvw2 


[5468 A] 


[4392 A] 


[3465 A] 


[2600 A] [2246 A] 


[1928 A] 


18.2 ± 0.1 


18.56 ± 0.07 


17.79 ± 0.06 


17.48 ± 0.06 17.50 ±0.06 


17.55 ± 0.05 



* Fixed value from measurements of the Galactic absorption by 
the Leiden- Argentine-Bonn survey (Kalberla et al. 2005). 

pile-up was found. Data were processed and screened by 
using the xrtpipeline task with default parameters and 
grades 0-12 (the single to quadruple pixels events). The 
spectrum was rebinned to have at least 30 counts per 
bin. 

No detection was found at E > 10 keV with BAT, after 
having binned, cleaned from hot pixels, deconvolved and 
integrated all the data available in this pointing. 

UVOT observed the source with all the six available 
filters. Data were integrated with the uvotimsum task 
and then analyzed by using the uvot source task, with a 
source region radius of 5" for the optical filters and 10" 
for the ultraviolet, while the background was extracted 
from a circular region l'-sized and centered in a nearby 
source-free region. It was not possible to select an an- 
nular region centered on PMN J0948+0022, because of 
nearby sources. The observed magnitudes were dered- 
dened according to the extinction laws of Cardelli et al. 
(1989) with A v = 0.277 and then converted into flux 



densities according to the standard formulae and zero- 
points (Poole et al. 2008). 
The results are summarized in Table 1. 

3.3. Radio 
3.3.1. Effelsberg 

The centimeter spectrum of PMN J0948+0022 was ob- 
served with the Effelsberg 100 m telescope on 24 January 
2009 (MJD 54855.5) within the framework of a Fermi re- 
lated monitoring program of potential 7-ray blazars (F- 
GAMMA project, Fuhrmann et al. 2007). The mea- 
surements were conducted with the secondary focus het- 
erodyne receivers at 2.64, 4.85, 8.35, 10.45, 14.60 and 
32.00 GHz. The observations were performed quasi- 
simultaneously with cross-scans, that is slewing over the 
source position, in azimuth and elevation direction, with 
adaptive numbers of sub-scans for reaching the desired 
sensitivity (for details, see Fuhrmann et al. 2008; An- 
gelakis et al. 2008). Pointing off-set correction, gain 
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Fig. 3. — The evolution of the radio spectrum of 
PMN J0948+022. Filled circles denote the Effelsberg obser- 
vations of January 2009. Archival data is shown in grey: triangles 
represent VLBA measurements conducted in October 2003 (Doi 
et al.2006), diamonds represent archival Effelsberg measurements 
obtained in 2006 (Vollmer ct al. 2008). 

correction, atmospheric opacity correction and sensitiv- 
ity correction have been applied to the data. 

The acquired radio spectrum has a convex shape with 
a turnover frequency between 10.45 and 14.60 GHz 
(Fig. 3). The low- frequency part spectral index o4°g| 5 , 
measured between 2.64 and 10.45 GHz, is -0.18 ± 0.02, 
whereas the high-frequency optically thin spectral index 
all 6 = 0.39. 

The comparision of the acquired spectrum with previ- 
ously observed ones (Vollmer et al. 2008, Doi et al. 2006) 
reveals that the source is presently in a much lower flux 
density state (see Fig. 3). This indicates intense vari- 
ability. From the change of the 4.85 GHz flux density 
over about three years (Vollmer et al. 2008) we estimate 
a variability brightness temperature (e.g Fuhrmann et 
al. 2008) of 1.4 x 10 11 K. Assuming the equipartition 
brightness temperature limit of ~ 10 11 K (Readhead et 
al. 1994), we obtain a lower limit for the Doppler factor 
of 5 > 1.4. 

3.3.2. Owens Valley Radio Observatory (OVRO) 

PMN J0948+0022 has been observed regularly from 4 
September 2007 at 16:25 UTC to 11 February 2009 06:53 
UTC (MJD 54347.68-54873.29) at 15 GHz by the Owens 
Valley Radio Observatory (OVRO) 40 m telescope as 
part of an ongoing Fermi blazar monitoring program of 
all 1159 CGRaBS blazars north of declination —20 de- 
grees (Healey et al. 2008). Flux densities were measured 
using azimuth double switching as described in Read- 
head et al. (1989). The relative uncertainties in flux 
density result from a 5 mJy typical thermal uncertainty 
in quadrature with a 1.6% systematic uncertainty. The 
absolute flux density scale is calibrated to about 5% us- 
ing the model for 3C 286 by Baars et al. (1977). This 
absolute uncertainty is not included in the plotted errors. 

PMN J0948+0022 has been reported to show variabil- 
ity by a factor of 2 in radio over year time-scales (Zhou et 
al., 2003) and ~ 31% fluctuations in month time-scales 
(Doi et al., 2006). The OVRO 40 m 15 GHz time series 
shows clear structure at timescales down to weeks, and 
year-scale fluctuations by a factor of 4 (Fig. 2, panel C). 
The rapid variability we observe in this object - at least 



400 mJy in 77 days or 5 mJy/day - enables us to deter- 
mine a variability brightness temperature of ~ 2 x 10 13 K, 
assuming the ACDM cosmology described in Sect. 1. 

It is often not easy to determine the optically thin spec- 
tral index of blazars at radio frequencies because they are 
complex structures, with different regions becoming op- 
tically thin at different radio frequencies. In the present 
case the most recent results show a turnover between 
10 and 15 GHz, and a 15 — 30 GHz spectral index of 
~ 0.4, but we do not believe that this is the optically 
thin spectral index. It is much more likely that one is 
still seeing synchrotron self-absorption so that the spec- 
trum between 15 GHz and 30 GHz is much flatter than 
the true optically thin spectral index. In such cases, it is 
safer to assume an optically thin spectral index of 0.75 
and to use the frequency of observation. These only have 
a small effect on the derived T cq unless a is very close to 
0.5, which is too flat, in our view, for an optically thin 
spectral index in most cases. 

The equipartition brightness temperature (Readhead, 
1994), in the current cosmological model, is then T eq ~ 
5.5 x 10 10 K, assuming an average optically thin spectral 
index of 0.75, and hence the equipartition Doppler fac- 
tor is 5 ~ 7, which is typical of highly variable blazars. 
This agrees with Doi et al. (2006), who reported an 
equipartition Doppler factor S > 2.7, and the Effelsberg 
lower limit reported in Sect. 3.3.1. This suggests that, 
in the compact radio emission regions in this object, the 
Lorentz T factor is of order 10. 

4. SPECTRAL ENERGY DISTRIBUTION (SED) 

Fig. 4 displays the spectral energy distribution (SED) 
built with the Fermi /LAT, Swift, and Effelsberg and 
OVRO data analyzed in the present work (red symbols) 
together with archival data (green symbols). Archival 
radio data are from Bennett et al. (1986), Becker et 
al. (1991), Gregory & Condon (1991), White & Becker 
(1992), Griffith et al. (1995), Doi et al. (2006); opti- 
cal/IR data are from USNO Bl for B, R, I filters (Monet 
et al. 2003) and from 2MASS for J, H, K filters (Cutri 
et al. 2003). 

The resulting SED strongly resembles that of a typical 
high power blazar, with two non-thermal emission peaks 
in the far IR and between 10 22 - 10 23 Hz (40-400 MeV). 
Also the peak produced by the accretion disk is well de- 
fined, due to the UVOT data, extending the photomet- 
ric coverage to the near UV. If we assume a standard 
Shakura & Sunyaev (1973) disk emission, the UVOT 
data permit fixing a lower limit to the mass of the black 
hole, which turns out to be around 1O 8 M0, in agreement 
with the estimates by Zhou et al. (2003). For lower 
masses, in fact, the luminosity needed to fit these data 
becomes super-Eddington. 

We model the contemporaneous optical to 7-ray data 
with a one-zone synchrotron and inverse Compton 
model, in accordance with the models generally used for 
blazars. Also for this source, the radio emission is as- 
sumed to come from larger-scale emission regions further 
away along the jet, while the rest of the SED is attributed 
to the region at the beginning of the jet plus the contri- 
bution of the accretion disk. The IR radiation produced 
by the assumed torus does not influnce the derived non- 
thermal SED, since the corresponding radiation energy 
density is much smaller than the one produced by the 
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Fig. 4.— Spectral Energy Distribution of PMN J0948+0022. Fermi/LAT (5-months data); Swift XRT and UVOT (5 December 2008); 
Effelsberg (24 January 2009) and OVRO (average in the 5-months of LAT data, indicated with a red diamond) are indicated with red 
symbols. Archival data are marked with green symbols. Radio data: from 1.4 to 15 GHz from Bennett et al. (1986), Becker et al. (1991), 
Gregory & Condon (1991), White & Becker (1992), Griffith et al. (1995), Doi et al. (2006). Optical/IR: USNO Bl, B, R, I filters (Monet 
et al. 2003); 2MASS J, H, K filters (Cutri et al. 2003). The dotted line indicates the contributions from the infrared torus, the accretion 
disk and the X-ray corona. The synchrotron (self-absorbed) is shown with a small dash line. The SSC and EC components are displayed 
with dashed and dot-dashed lines, respectively. The continuous line indicates the sum of all the contributions. 



lines. 

The complete description of the general model used 
can be found in Ghisellini & Tavecchio (2009); here we 
briefly summarize the main parameters: 

• the source is assumed to be located at a distance 
-Rdiss = 6.75 x 10 16 cm from a black hole of mass 
M = 1.5 x 1O 8 M ; 

• the source is assumed to be a sphere, of radius R = 
6.75 x 10 15 cm (i.e. a conical jet of semiaperture 
angle tp = 0.1 rad is assumed); 

• it moves with a bulk Lorentz factor T = 10 at a 
viewing angle 8 V = 6°; 

• the disk luminosity is 40% of the Eddington value 
(L disk = 9 x 10 45 erg s" 1 ); 



• above (and below) the disk we assume an X-ray 
emitting corona, producing a luminosity L cor = 
0.3Ldisk and with a power law spectrum (energy 
index ax — 1), ending in an exponential cut at 
150 keV; 

• 10% of the disk emission is assumed to be absorbed 
and re-emitted by the broad-emission line region 
(BLR), which emits narrow permitted lines, in this 
specific case; 

• the BLR is placed at a distance -Rblr = 
10 17 L^( s 2 k 45 = 3 x 10 17 cm, where Ldisk,45 is the 
disk luminosity in units of 10 45 erg s _1 ; 

• a dusty torus absorbs 10% of £disk re-emitting 
it in the far-IR. Its distance is assumed to be 
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Fig. 5. — Radio versus 7-ray properties of PMN J0948+0022 (indicated with a filled square) compared with the other blazars detected 
by LAT (figure adapted from Fig. 14 in Abdo et al. 2009b). (left panel) Peak 7-ray flux (E > 100 MeV) vs. radio flux density at 8.4 GHz; 
the dashed lines show the CRATES flux density limit and the typical LAT detection threshold, (right panel) 7-ray photon index vs. radio 
luminosity. 



R 



in. 



2.5 x 10 18 Lj{ s 2 k 45 cm, calculated by con- 



sidering that the dust is at the temperature of 370 
K (see Cleary et al. 2007), where is the maximum 
of efficiency in reprocessing the impinging photons 
into IR radiation. 

For the above choice of parameters, the BLR is the 
main producer of the seed photons for the inverse Comp- 
ton scattering in the emitting region. In Fig. 4, this 
component is labeled EC. On the other hand, the Syn- 
chrotron Self-Compton (SSC) radiation is important in 
the X-ray band, where it dominates the radiative output 
(the two-peak shape of the SSC curve correspond to the 
first and second order scattering). In this scenario, i.e. 
with the X-ray data due to the SSC emission, the mag- 
netic field B is strongly constrained, since it controls the 
value of the SSC emission. In our case B = 3.2 G. 

The particle energy distribution of the emitting elec- 
trons is calculated by the continuity equation, assuming 
the injected electrons are distributed in energy according 
to a smoothly joining broken power law of slopes 7 _1 and 
7~ 2 - 2 below and above 7 = 800, respectively. The max- 
imum electron energy corresponds to a random Lorentz 
factor of 7 m ax = 1600 and the inverse-Compton peak cor- 
responds to 7 poa k = 411. The total power injected in the 
form of relativistic electrons is L' c = 2.4 x 10 43 erg s _1 , 
as calculated in the comoving frame. 

Once accelerated, the bulk kinetic power carried by 
the electrons in the jet is L e = 5.0 x 10 44 erg s _1 , to be 
compared with the Poynting flux Lb = 1-8 x 10 44 erg s _1 
and the power in radiation, L ra d = 2.0 x 10 45 erg s _1 . 
As occurs in typical powerful blazars (see e.g. Celotti 
& Ghisellini 2008, also for the exact definitions of these 
powers), the radiation observed carries more power than 
what is available in magnetic field and electron energy 
Thus we require that the jet also carries protons: if we 
assume one proton per electron, we obtain L p = 4.8 x 



the presence of pairs. 

The impact of the presence of narrow permitted lines 
in the BLR can be evaluated in two ways. First, there 
can be a geometrical explanation: if the BLR is torus- 
like and we are observing it face-on, then the Doppler 
broadening is reduced and the line width is smaller than 
usual. A similar geometry has been invoked by Decarli et 
al. (2008) to explain why NLS1 have not small masses. 

This BLR geometry, although different from what used 
in the adopted model, does not affect any of the model 
results, since the important parameter is the angle at 
which the jet sees the BLR. Second, it is still possible to 
consider a isotropic geometry of the BLR and take into 
account the effect of the radiation pressure, as proposed 
by Marconi et al. (2008). The BLR is then moved farther 
off the central singularity (-Rblr = 5 x 10 cm), but 
the impact in the model parameters mainly results in 
an increase of F from 10 to 13. Other new values of 
parameters are: 7 max = 1500, 7 pC ak = 430 L' c = 2.6 x 
10 43 erg s _1 , B — 2.01 Gauss. The power output of the 
jet is: L c 



1.2 x 10 



45 



erg s 
1 



Lt 



1.3 x 10 44 erg s"\ 



10 46 



erg s . This value can be reduced, by assuming 



i^d = 3.2 x 10 45 erg s" 1 , L p = 7.9 x 10 46 erg s" 1 . 

It is worth noting that the systematics in the LAT data 
do not affect significantly the results of the modeling. 
Indeed, in the case of a LAT flux equal to the 130% of 
the observed value (the worst case), it is sufficient to 
increase a little the injected power L' e (from 2.4 x 10 43 to 
2.8 x 10 43 erg s _1 ) and to decrease a little the magnetic 
field B (from 3.2 to 2.5 Gauss). The derived quantities 
become L e = 6.3 x 10 44 erg s^ 1 , Lb = 1.0 x 10 44 erg s _1 , 
L rad = 2.5 x 10 45 erg s" 1 , L p = 5.0 x 10 46 erg s" 1 . A 
change in the 7-ray photon index equal to ±0.1 results 
in very negligible changes in the SED. 

The model parameters used to explain the SED of 
PMN J0948+0022 are rather well constrained, given the 
basic assumptions of the model. However, they may not 
be unique. Another solution might be possible if we as- 
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sume that the dusty torus contributing to the IR seed 
photons is hotter and more compact than assumed here. 
It must be hotter because otherwise we need too ener- 
getic electrons to fit the hard X-ray emission (by EC), 
and these very same electrons would overproduce (by 
synchrotron) the observed emission in the optical. A 
stratified and clumpy torus (e.g. Ncnkova ct al. 2008) 
could work. Also in this "hot and clumpy torus" case 
we can obtain a reasonable representation of the data 
(i.e. the same ratios between the radiation energy den- 
sities - both by synchrotron and external photons - and 
the magnetic energy density), at distances one order of 
magnitude larger than assumed here. As a consequence, 
the emitting region should be one order of magnitude 
larger, and thus should vary on longer timescales. The 
detection of a typical variability timescale can therefore 
discriminate between the two possible solutions. 

5. DISCUSSION AND CONCLUSIONS 

Our findings show clearly that PMN J0948+0022, the 
first narrow-line quasar detected in 7-rays, hosts a rela- 
tivistic jet with SED very similar to those of "classical" 
FSRQ. 

The parameters relevant for the spectral shape of the 
SED derived from the modelling, using the approach by 
Ghiscllini and Tavecchio (2009), are consistent with those 
of high-power blazars, though at the lower limit in power 
of the FSRQ region as derived, for example, in Celotti & 
Ghiscllini (2008) for a large sample of blazars (cf. their 
Fig. 2 and Fig. 6). The comparison of radio vs 7-ray 
properties shows that in both bands PMN J0948+0022 
has a relatively low power with respect to the other 
FSRQ detected by LAT to date, although no striking 
differences appear (Fig. 5). Effelsberg and OVRO radio 
observations also show characteristics similar to those of 
FSRQ. 

We believe that PMN J0948+0022 could be one of the 
first examples of a FSRQ with relatively small mass but 
high accretion rate in terms of Eddington ratio. This 
would fit well with the scenario in which the transition 
between high peaked and low peaked blazars is not re- 
lated to absolute power but to the Eddington ratio as 
indicated by the FRI-FRII separation which depends on 
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